BOUNDS ON ORDERED CODES AND ORTHOGONAL ARRAYS 



ALEXANDER BARG* AND PUNARBASU PURKAYASTHAt 



Abstract. We derive new estimates of the size of codes and orthogonal arrays in the 
ordered Hamming space (the Niederreiter-Rosenbloom-Tsfasman space). We also show 
that the eigenvalues of the ordered Hamming scheme, the association scheme that describes 
the combinatorics of the space, are given by the multivariate Krawtchouk polynomials, and 
establish some of their properties. 
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1-A. The Niederreiter-Rosenbloom-Tsfasman metric space. Let Q be a finite alphabet 
of size q viewed as an additive group mod q. Consider the set Q'"'" of vectors of dimension 
rn over Q. A vector x will be written as a concatenation of n blocks of length r each, 
X = {xii, . . . , xir', . . . ; Xni, . . . , Xnr}- For a given vector x let ei,i — 1, . . . , r be the 
number of r-blocks of x whose rightmost nonzero entry is in the ith position counting 
from the beginning of the block. The r-vector e — (ei, . . . , e^) will be called the shape 
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of X. For two vectors x,y E Q'"'" let us write a; ~e 2/ if shape(a: — y) = e. A shape 
vector e ~ (ei, . . . , Cr) defines a partition of a number N < n into a sum of r parts. Let 
Ar,n = {e G U {0})'' : J2i — be the set of all such partitions. For brevity we 
write 



Let X £ Q'' " be a vector of shape e. Define a weight function (norm) on Q'' " by 
setting w{x) — |e|' and let dr{x, y) = w{x — y) denote the metric induced by this norm. 
We call the function dr the ordered weight. The ordered weight was first introduced by 
Niederreiter |[29 | in his study of low-discrepancy point sets. Later, Rosenbloom and Tsfas- 
man f34l independently defined the weight w(a;) (more precisely, the weight w defined 
below), calling it the m-metric, and studied codes in Q'''" with respect to it. The set Q'''" 
together with this metric will be called the ordered Hamming space (the NRT space) and 
denoted hy H — H{q, n, r). Unless specified otherwise, below by distance (weight) we 
mean the ordered distance (weight) for some fixed value of r. Note that the case r = 1 
corresponds to the usual Hamming distance on Q". 

An {nr, M, d) ordered code C C Q'^'" is an arbitrary subset of M vectors in Q'"'" such 
that the ordered distance between any two distinct vectors in C is d or more. If Q is a finite 
field and C is a linear code of dimension k, we refer to it as an [m\ k, d] code. 

To define ordered orthogonal arrays, let us call a subset of coordinates I C {!,..., rn} 
left-adjusted if with any coordinate ir + j,0 < i < n — 1,1 < j < r in the ith block 
it also contains all the coordinates {ir + 1, ir + j — 1) of the same block. A subset 
A C Q*"'", \ A\ — M is called a [t, n, r, q) ordered orthogonal array (OOA) of strength t 
if its projection on any left-adjusted set of t coordinates contains all the rows an equal 
number, say 9, of times. The parameter 9 is called the index of A. It follows that M = 9q*^ . 
If C is a linear [nr, k , d] code then the code orthogonal to it with respect to the dot product is 
a (d — 1, n, r, q) linear OOA of index 9 = jf Q is equipped with the structure 

of an additive group, then one can construct additive OOAs. If Q is a finite field, it is 
possible to construct linear OOAs. 

OOAs (also called hypercubic designs) were introduced in Lawrence ll20l and Mullen 
and Schmid 128 i i as a combinatorial equivalent of point sets suitable for numerical inte- 
gration over the cube. Informally this link can be described as follows. Let C„ = [0, 1]" 
be a unit cube, let / be a continuous function of bounded variation and let A/" be a set of 
M ~ points in C„, called a net. It is known that the error of quasi-Monte Carlo inte- 
gration I fdx — g^™ X^xGA/' b^ bounded above by V{f)D{J\f) where V{f) 
is the total variation of / on C„ and D{Af) is the discrepancy factor of the net. The pa- 
rameter D{M) measures the deviation of the net from a uniformly distributed set of points. 
The study of uniformly distributed point sets was initiated by H. Weyl and E. Hlawka. 
Sobol |36J developed the notions described above and gave the first constructions of nets 
with bounded factor D. The study of nearly uniform point sets was taken up in Niederre- 
iter |291 which put forward the notion of (i, m, s)-nets and derived a bound on D{Af) via 
its parameters. We refer to Il23l[3ni30l for detailed background and more references for 
{t, m, s)-nets and in particular, to the literature on their constructions. 

The following theorem relates OOAs and {t, m, s)-nets. Below we use its statement as 
a definition of a (t, ni, s)-net. 

Theorem 1.1. (Lawrence |20i , Mullen and Schmid ||28l). There exists a (t, m, s)-net if 
and only if there exists an (m — t, s,m — t, q) OOA of index g* and size M = g™. 
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Independently of this line of work, Rosenbloom and Tsfasman ll34l considered codes 
in the ordered Hamming space and derived bounds on their size. It became clear shortly 
after their work that OOAs and ordered codes are dual types of objects in the sense of 
Delsarte's algebraic theory of coding. This link opened up an avenue for applications of 
coding-theoretic methods to the study of {t, m, s)-nets and motivated the study of ordered 
codes and OOAs independently of these applications. In particular, Martin Il24ll25ll22ll26l 
has constructed an association scheme that describes the combinatorics of the space H, 
formulated the linear programming bound (LP bound), outlined a construction of orthogo- 
nal polynomials that describe the eigenvalues of the scheme and derived Plotkin and Rao 
bounds on OOAs via linear programming. Much attention was also devoted to relations 
between the weight enumerators of linear codes and their duals in the ordered Hamming 
space and its generalizations. In particular, a Mac Williams theorem for the NRT space was 
derived by Martin |i24|, see also Skriganov |35|, Dougherty and Skriganov [ 14|. Bierbrauer 
ESUni studied coding constructions of {t, m, s)-nets and bounds for them including the 
LP bound. Apart from the combinatorial context, ordered codes arise in a number of ap- 
plied problems such as recent algebraic list decoding algorithms of Reed-Solomon codes 
II32I . a study of linear complexity of sequences lIZTl . and in a problem in communication 
theory 

In this paper, we derive several new bounds on OOAs and ordered codes. We begin 
with a Bassalygo-Elias bound on codes which improves the known upper bounds on their 
size. However, the bulk of our results are devoted to code bounds using the approach via 
association schemes and linear programming. We begin with a study of the eigenvalues of 
the ordered Hamming scheme, an association scheme that describes the combinatorics of 
the space H. The approach that we follow relies on the orthogonality relation of the eigen- 
values. This enables us to identify the eigenvalues as multivariate generalizations of the 
well-known Krawtchouk polynomials, i.e., a family of real polynomials of r discrete vari- 
ables orthogonal on Ar,„ with respect to the weight given by the multinomial probability 
distribution. This is the subject of Section|3] Turning to bounds, in Sect.|5]we derive a new 
universal estimate of the size of ordered codes and OOAs with a given distance (strength). 
The asymptotic version of this estimate improves the other results in a certain range of 
rates. The final section is devoted to the case r = 2 for which the bounds can be further 
improved relying on a du^ect approach. 

1-B. Notation. Together with the space H we will consider the space H{q, n, r) which 
differs from it in that the vectors are read "from right to left." Namely, for x £ Q^'^ 
let shapc(a;) = (ei , . . . , e^), where Cj is the number of blocks {xn , ■ ■ ■ , xtr) such that 
Xii — ■ ■ ■ — Xi,r-j — and Xi,r-j+i 7^ 0. Let w(a;) = |e|', where e = shape(a;), and let 
dr{x,y) = w(a; - y). 

These metric spaces are identical to each other; the reason for considering them both 

— > 

is that if we equip H with the structure of an association scheme then its dual scheme in 
Delsarte's sense gives rise to the space H . In particular, if Q is a finite field and C is a 
linear code in H then its dual code lives in H . We elaborate on this below. _^ 

An easy combinatorial calculation shows that the number of vectors of shape e in is 
given by 



(1) 
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and the number of vectors of weight d equals to 

(2) Sd= ^ Ve- 

e: \e\'—d 

Let A{z) = (q — l)z{z^ — l)/{q{z — 1)) and let zq = zq(x) satisfy the equation 

— 1 ^ 

(3) xr(l + A(z))^- y"iz\ 

Define the function 

Hq^r{x) = x{l - \ogg zo) + ^ log,(l + A{zn)) ■ 
In the case r = 1 we write hq{x) instead of Hq,i{x), where 

X 

hg{x) ^~x\og ---(l-a;)log 

^ g — 1 ^ 

Let 

1 q'' -1 



(4) 6„n^l--J2l 



1 



r rq^ <? — 1 

The asymptotic volume of the sphere in H is given in the next lemma. 

Lemma 1.2. II34II (a) For < a; < 1, equation Q has a unique positive root zq{x), zq G 
[0,r]. 

(b) Let r > 1 be fixed, n — > oo, d/n — > r5, then 

<5crit <S <1. 



(5) lim (nr)-i log, V 5, = <^ 



1=0 

2. Bounds on Ordered Codes and Arrays 

In general, given the value of the distance or the strength, our goal is to construct as 
large codes and as small OOAs as possible. The latter will also account for small-size 
(t, m, s)-nets with bounded discrepancy. In this section we recall the known bounds on 
ordered codes and OOAs and derive a new bound on the size of codes. 



2- A. Existence bounds. A Gilbert-type bound on ordered codes was derived in 113411 . 
Theorem 2.1. There exists an {nr, M, d) code in the space H whose parameters satisfy 



1=0 

If Q is a finite field, then there exists a linear code with the same parameters. 

A bound that applies specifically to linear codes was proved in fSl. It is analogous to 
the Varshamov bound for the Hamming space. 

Theorem 2.2. Suppose that m and t satisfy the conditions 

t-T 

^^,,„_i <g"-(^-i), r=l,...,t-l. 

Then there exists an [nr, nr — m] linear code in H of distance > i + 1, and a (t, n, r, q) 
linear OOA of dimension m. 
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2-B. Nonexistence bounds. While in general bounds on codes do not imply lower bounds 
on OOAs, there are two special cases when these two types of results are equivalent. First, 
if C is an [nr, fc, d\ linear code over F^*" then the code C :~ {y £ F^*" : X)"=i ^iVi ~ 
for all X E C} is a (d — 1, n, r, q) linear OOA. Next, if an upper (lower) bound on codes 
(OOAs) is obtained by linear programming as explained in the next section then the same 
solution of the LP problem gives a lower (upper) bound on OOAs (codes). 
We next mention some upper bounds on codes and OOAs. 

2-B.l. Singleton bound. The parameters of any {nr, M, d) code satisfy 



2-B. 2. Plotkin bound. A Plotkin bound on codes was established in 113411 . Namely, the 
following result holds true. 

Theorem 2.3. Let C <Z H be a code of size M and distance d > nrScni- Then 

d 



M < 



d - nrScrii 

A dual Plotkin bound on OOAs was proved by Martin and Visentin 
Theorem 2.4. ||26l Let C be a [t, n, r, q] OOA. Ift > nrSait - 1 then 

nrScrit 
t + 1 



\C\>q"^ll- 



2-B. 3. Hamming-Rao bound. According to the Hamming bound, the parameters of any 
{nr, M, d — 2t + 1) code satisfy 



rn 

M< ^ 



A dual bound in this case is the Rao bound which for the NRT space was established by 
Martin and Stinson ||25l : The size M of any {t — 2t, n, r, q) OOA satisfies 

r 
1=0 

2-B. 4. A Bassalygo-Elias bound on codes. The next result is new. 
Theorem 2.5. Let C be an {nr, M, d) code. Then 

1 



M < q''"dn 



Sw{dn - 2wn + 



for any w < nr6cnt{l - ^1 - d/ {nr dcdt))- 
Proof: We will rely upon the next lemma. 

Lemma 2.6. Let C C H , \C\ ^ M be a code all of whose vectors have weight w and are 
at least distance d apart. Then for d > 2w — w'^ / {nrS^nt), 

dn 



M < 



dn - 2wn + 
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Proof: Let C* be a projection of C on the ith block of coordinates. For a vector z E Q'' let 
z'^ = (zr-h+i, . . . , Zr) be its suffixof length ft,. Given a; e C, we denote by a;* G its ith 
block and write a;*''' to refer to the ft-suffix of x\ For i = 1, . . . , n; /i = 1, . . . , r; c e Q'' 
let 

be the number of vectors in the zth block whose /^-suffix equals c. We have 

r 
r 

Compute the sum of all distances in the code as follows: 

n r 

dr{x,y)^nrM'-j2 E E E 

x,yeC i=l x\y^£C^ h=l cGQ'' 



(6) =nrM'-Y^Y.T.i^t)"- 

1=1 h=i ceQi^ 

To bound above the right-hand side, we need to find the minimum of the quadratic form 

n r n r 

i=l /i=lcee''\{0} i=l h=l 

under the constraints 

n r 

EE ^^.0 = M{nr - w), ^ = A/ (1 < < r, 1 < z < n). 

i=l h=l ceQ*^ 

Critical points of F in the intersection of these hyperplanes, together with (|7]), satisfy the 
equations 



1 < i < n; 1 < /i < r;c e Q''\{0} 



(8) 2^t + A.-0, 

2X'}jy + a + I3i^h = 0, I <i <n;l <h <r 



The system (I?])-® has a unique solution for the variables ^, j3i,in ol\ in particular, 



1 



9 

Mw 



- 1 



i.c h ' 

q'^nr d„ii 



/i = 1, . . . ,r,z = 1, . . . ,n, 



/i = l,...,r,z = l,...,n,ce Q''\{0}. 



To verify that this critical point is in fact a minimum, observe that the form F is convex 
because its Hessian matrix is 2/ and is positive definite (both globally and restricted to the 
intersection of the hyperplanes (|7]i ). Substituting these values of the As and taking account 
of the fact that q^'' = r(l — (5crit), we get 



n r 

^^EEEt 

i=l h=l c^O 



Mw \2 



q^nrSc 



i h 

<? 2w 



w 



nrSn, 



+ 1 



cnt 



n 
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Then from (|6| we obtain 

M2 



rfM(M-l)< dr{x,y) < —(^ 



2wn — 



n \ r(5c, 

x,y£C 



which gives the result. 



The proof of the theorem is completed as follows. Let C Q*^'" be a sphere of radius 
w around zero. Clearly, 

\C\\S^\ = J2 \{C-x)nS^\<q''^Ag{nr,d,w), 

where Aq{nr, d, w) is the maximum size of a distance-d code in Sw With the previous 
lemma, this finishes the proof. I 

Remarks. 1. This theorem implies a lower bound on the size M of a linear OOA 

{t - 1, n, r, q): for any w < nr4rit(l - ^1 - t/{nrScrit)), 

(9) M >—S^ (tn - 2wn + 



tn V rJci-it 

and in particular, a lower bound on linear (m — r, to, n)-nets, m = log^ M. 

2. Caution should be exercised in dealing with codes of a constant weight in the NRT 
space, i.e., codes on the sphere S^ va H . Indeed, the sphere Si^ together with the metric 

is not homogeneous: in particular, the number of points in located at a given distance 
from a point x & Sw depends on x. However, this does not cause problems in the previous 
theorem. 

3. The argument used in the proof of Lemma l2!6] can be also used to give a proof of the 
Plotkin bound, Theorem l2.3l that is simpler than the ones known in the literature. Indeed, 
let C C iJ be a distance-d code. Consider again expression (|6]l and note that this time there 
is no restriction on the weight of the codewords. Using the Cauchy-Schwarz inequality and 
the fact that X^ceF'' -^i'c — obtain 

M{M - l)d < nrM^ " XI XI ~ = nrM'^Scnt- 

1=1 h=l ^ 

Solving for M concludes the proof. 

2-B.5. Asymptotics. In this section we assume that n — > oo and r is a constant. For a code 
of size M let i? = log^ M/nr be the code rate. Given a sequence of (rrii, Mi, di) codes 
we will say that its asymptotic rate is R and the asymptotic relative distance is 5 if 

lim ^— log Mi — R, lim = S. 

i^oo rui ^ i-»oo rui 

The Plotkin bound implies that the asymptotic rate and distance of any sequence of codes 
satisfy 

c 

Ocrit 



i? = 0, s>s, 



cnt- 



g 
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To state the "sphere packing" or "volume" bounds on ordered codes we rely upon 
Lemma [TT2I Namely lf34l . there exists a sequence of [rni,ki,di] linear codes Ci,i = 
1,2,..., such that rii — > 00, ki/{rni) — > R, di/{rni) — > 5 and 

i? > 1 — Hq r{S), < 6 < Sa-n (Gilbert- Varshamov bound) . 

On the other hand, for any such sequence of codes, 

i? + iJq.rCV^) < 1 (Hamming bound) . 

The asymptotic version of Theorem l2.5l is as follows: 

Theorem 2.7. (Asymptotic Bassalygo-Elias bound). For < 6 < 6cnt the asymptotic rate 
and distance of any sequence of codes satisfy 

(10) R<1- Hq.r (4rit(l - - 5 /Sent)) ■ 

This bound is better than the Hamming bound for all 5 e (0, Jcrit] • It is also often better 
than the Plotkin bound. For instance, for q = 2,r = 2 the bound (fTOl i is better than the 
Plotkin bound for all 5 G (0, 5„ii)- For larger g, r the improvement is attained only for low 
values of 5 since the right-hand side of (flOl i becomes n-convex close to S^^^. For instance, 
for g = 3, r = 4 this range is (0, 0.54), etc. 

2-B.6. Asymptotic bounds for digital (t, m, s)-nets. A it, m, s)-net is called digital if the 
OOA that corresponds to it forms a linear subspace of F^*". Therefore, bounds on linear- 
OOAs apply to the special case of digital (i, m, s)-nets. However, studying asymptotics 
for this case requires a different normalization since the strength ra — t of the OOA that 
corresponds to the net equals r, and both approach infinity independently of s. Therefore, 
let i? = m/ s denote the rate and 5 = {m—t)/s denote the relative strength of the OOA that 
corresponds to the net. To state the bounds, we need to compute the asymptotic behavior of 
the volume of the sphere, which is different from (|5]l. The next result is due to Bierbrauer 
and Schmid ID. 

Theorem 2.8. There exist families of digital {t, m, s)-nets with s, (m — t)^ 00 for which 
(R, 6) satisfy the bound R < 5* ((5), where 

and a is defined by Sa{q — 1 + a) = (q — — a). 

On the other hand, by the Rao bound, any family of {t, m, s)-nets satisfies R> '^{S/2). 
Observe that Theorem l2.5l in this case gives the same result as the Rao bound because the 
increase of the packing radius in (|9]l over S/2 vanishes asymptotically. Indeed, taking 
uj = w/n and replacing t with m ~ t,we obtain from (|9]l 

M> l{S~2uj + o{l))S^n. 


The tightest bound is obtained if we take oj = S/2 in this inequality. 

Remark: We note that in the case that both n ^ 00 and r ^ 00 while S = d/nr tends 
to a constant bounded away from and 1, the lower and upper bounds on codes coincide 
|[34ll (the Gilbert- Varshamov bound converges to the Singleton bound). 
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3. Association Schemes, and the Ordered Hamming Scheme 

The coding-theoretic notion of duality prompted a study of the association scheme that 
describes the combinatorics of the NRT space. We briefly recall some elements of Del- 
sarte's theory of association schemes 1 13 1. A symmetric association scheme with D classes 
is a finite set X, |X| = N, equipped with a set TZ = {Rq, Ri, . . . , Rd} of symmetric bi- 
nary relations on X x X such that 

(i) Ro = {{x, x)}, R, n Rj = 0, ugoR, =XxX; 

(ii) For each < i, j, k < D the number 

pIj = \{ix, y) e Ri, {x, z) e Rj given that (y, z) e Rk}\ 

depends only on k). Moreover, p^^ — i- 

The parameters p^ ^ are called the intersection numbers of the scheme A — {X,Tl). 
The numbers Vi = Pii are called the valencies of A. Let Ai be the adjacency matrix of the 
relation Ri,i — 0, . . . , DAtis clear that 

D 

Aq = I, Ai — J (all-one matrix), 

i=0 

and for all < i,j < D the product AiAj is contained in the linear span of {Aq, Ai, 
. . . ,^d}- The matrices Ai form a commutative algebra over C with respect to matrix 
multiplication, called the Bose-Mesner algebra. 

It is clear that the Bose-Mesner algebra is also closed under the Hadamard (element- 
wise) multiplication o, viz. Ai o Aj = SijAj. With respect to o, this algebra has a basis of 
primitive idempotents {Eq, Ei, . . . , E^} that satisfy 

1 1 ^ 

£;o = - J, E,oEj^-Y. '^t-^fc' 0<h3<D. 

The numbers qfj are nonnegative. They are called the Krein parameters of the association 
scheme. The quantities Hi = qfi are called multiplicities of the scheme A. The matrices 
P and Q defined by 

D 

and 

1 ^ 

^ j^Y^Q^jA^, 0<j<D, 

1=0 

are called the first and second eigenvalues of A, respectively. The eigenvalue matrices 
satisfy PQ = NI. Further, for < i, j < D the eigenvalues satisfy 

(11) HiPij = VjQji, 



(12) 



(13) 



^HkPkiPkj = NviSij, 

k=0 

D 



1=0 
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Two D-class association schemes are called Delsarte duals of each other if the adjacency 
matrices Ai, the first eigenvalues P, and the intersection numbers p'f ^ of one scheme are, 
respectively, the idempotents Ei, the second eigenvalues Q and the Krein numbers q^j of 
the other. The duality also exchanges the role of the matrix and Hadamard multiplication. 
A scheme is called self-dual if it equals its dual. For instance, the Hamming scheme 
Hn = {X ^ Q", Ri = {{x, y) e X, dnix, y) = i},i = 0, ... ,n) is self-dual. One of the 
manifestations of self-duality in this case, obvious from the definition, is the Mac Williams 
theorem that relates the weight distribution of an additive code to that of its dual. A scheme 
is called formally self-dual 1 11 1 if there exists some ordering of primitive idempotents 
under which P = Q. In a formally self-dual scheme Vi ~ fii and p'Hj — q^j. Following 
Delsarte lfT3l p. 17] a scheme {X",TZ) is called an extension of an r-class scheme K, — 
{X, T) = {Dq, Di, . . . , Dr)) if its vertex set is the 7i-fold Cartesian product of X and the 
relations i?e, e G A^.n are given by 

Re = {{{Xll, . . . ,Xi„), {X21, . . . ,X2n)) ■ 

\{j : {xij,X2j) e All =ei, i = 0,l,.. .,r}. 

Apart from fTsl we refer to lISl fTTlfTsll for the proofs of these results and more information 
on association schemes. 

The association scheme for the NRT space was constructed and studied by Martin and 
Stinson 1241 . Define an r-class "kernel scheme" IC{Q^'^,'D = {Dq, Di, . . . , Dr)) with the 
relations 

A = {{xi,X2) G Q"^i X Q'-'i : dr{xi,X2) = i}, i^O,l,...,r. 

Theorem 3.1. Il24l The space X = Q'"'" together with the relations 

Re = {{x, y) e X X X : X ^eV} (e G A^,,,) 
— ^ 

forms a formally self-dual association scheme Ti, called the r -Hamming scheme. It can be 
constructed as an n-fold Delsarte extension oflC. 

The dual scheme of TL is Ti. whose point set is X — Q*"'" and the set of relations is 
given by 

Re = {ix,y) e X X X : shape(a; - y) = e} (e G Ar,„). 
4. Multivariate Krawtchouk Polynomials 

In the conventional case of r = 1, eigenvalues of the Hamming scheme are given by the 
Krawtchouk polynomials 

(14) ^^(-'-)=E(-i)'('^-i)'"0(n) 

which form a family of polynomials of one discrete variable orthogonal on the set {0, 1 , . . . , ti} 
with weight a{i) — (")(q — i.e., the binomial probability distribution. Here we 

are interested in their generalization for the r-Hamming scheme. 

Observe that the valencies Ve = Pe e of scheme are given by ([T]i. By self-duality 
and (fTZt . the eigenvalues are orthogonal on the space of partitions Ar,„ with weight Ve. 
Below it will be convenient to normalize the weight. Let V — Vr^n be the space of real 
polynomials of r discrete variables x = (xi, X2, . . . , x^) defined on Ar.„. Let us define a 
bilinear form acting on the space V by 

(15) (mi,U2)= ^ Mi(e)M2(e)w(e), 
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where w{e) — q "''we- Letting pi — ^{q — — 1, . . . ,r; po = q '', we observe 

that 



■1=0 

forms a multinomial probability distribution on Ar,„. Therefore, r-variate polynomials 
orthogonal with respect to this weight form a particular case of multivariate Krawtchouk 
polynomials. 

For a partition / e ^r,n denote by 

Kfix) ^Kf,^,„j,.{xi,...,Xr) 

the Krawtchouk polynomial that corresponds to it. Let k = |/| be the degree of Kf. 

Our goal in this section is to derive properties of the polynomials Kf. In their large 
part, these properties are obtained by specializing to the current case general relations of 
the previous section. However, some work is needed to transform them to a concrete form 
which will be used in later calculations. 

The following relations are useful below. 

Lemma 4.1. 

(16) {x^,l) =niq^l)q'-^-\ z = l,...,r 

(17) {x^,XJ) ^n{n~l){q-lfq'+^-^'-^, l<t^J<r 

(18) {x^,x^) =n{q-l)q^—\l + {n-l){q-l)q'-'^-^), i = l,...,r. 
Proof: To prove ( fT6] l, compute 

e j — 1 

\tiO, cl, ■ ■ ■ ,Ci i, ■ • ■ , Cr/ 

^ J — 

The sum on e on the last line equals {q — which finishes the proof. The 

remaining two identities are proved in a similar way. I 

4- A. Properties of the polynomials Ke- 

(i) Ke{x) is a polynomial in the variables xi, . . . ,Xr of degree k = \e\. There are 

C^r^-T^) different polynomials of the same degree, each corresponding to a partition of k. 

(ii) (Orthogonality) Equation (fT2] l is rewritten as 

(19) {Kf,Kg)=VfSf^g, 11^/11 =v^- 

In particular, let Fi — (0'~^10''^*^^), i = 1, . . . , r be a partition with one part. We have 

(20) \\Kf£ ^ {Kf,,Kf,) ^n{q-l)q'~\ i^l,...,r. 

Indeed, equality ( fT9] l is simply ( fT2] i specialized to the case at hand and (|20] | is obtained 
from ([TJ. 

(iii) (Linear polynomials) For i — 1, . . . ,r, 

(21) Kpiix) = q'^^{q - l){n - Xr Xr-i+2) - q^Xr-i+i- 
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Proof: This is shown by orthogonalizing the set of linear polynomials (1, xi, X2, ■ ■ ■ , 
Xr). We take Kq^... ,o — 1- Use Lemma l4n to compute 

Kp^ix) = ci{xr - {xr, 1)) = Ci[xr - n{q - l)/q) 

for some constant ci. To find ci, use ( |20] |: 

n(g - 1) = c?||x. - ^ - l)^-^. 

Hence ci = zLq. We take Kpiix) ~ n{q — l) — qxr choosing ci = —qsothatKp-^{0) > 0. 
Next let us perform the induction step to compute Kpt^i (x): 

i 

(22) Kf,^,{x) - C,+,(xr-^ - \\KF,r'{Xr-^,KF^)KF^{x)y 

3=0 

where the polynomials Kf , j ^ 0, . . . ,i, have the form (ISTT l by the induction hypothesis. 
Straightforward calculations using (fT6b-(fT8b show that 

^Fi+iix) = Ci+i{Xr^i - {{q - l)/q){n - Xr Xr-i+l))- 

Again using (|20|) . we find that q+i — ±(7'+^ ; as above, we will choose the minus. I 

(iv) The next property is a special case of (fTTT l. 

VeKf{e)=VfK,if) (e,/eA,,„). 

In particular, 

(23) Kf{0) = vf. 

(v) For any f,g e A^,™ 

(24) Kf{e)Kg{e)= ^ 

where the linearization coefficients g = \{z E Q''^" : z ^ f x, z y;x y}\ are 
the intersection numbers of the scheme. In particular, ^ > 0. This is a special case of 
property (fT3] l. 

(vi) (Three-term relation) Let IK„ be a column vector of the polynomials Kf ordered 
lexicographically with respect to all / that satisfy |/| = k. The three-term relation is 
obtained by expanding a product P{e)Kii{e) in the basis {Kj}, where P{e) is a first- 
degree polynomial. By orthogonality, the only nonzero terms in this expansion will be 
polynomials of degrees k + 1, k, k — 1 ifTSl p. 75]. 

We will establish an explicit form of the three-term relation for P(e) = (5ciit?'fi — |e|'. 
We have 

(25) P{e)K^{e) - a,K,+i{e) + 6«K,(e) + c,IK,-i(e), 

where a^, b^, c„ are matrices of order {^^^^^) x (J'^^-i^^) '^^^ * = 1, 0, — 1 , respectively. 
The nonzero elements of these matrices have the following form: 

a^[f,h]=L,{n + l) if /i= (/i,..., /, + !,...,/,), 

c4f,h]^L,{n-K+l)q'-^{q-l) if = (/i, ...,/,- 1, 
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(26) 

where L; — 



{ UM-\q-2) if h = f, 

Hfk + l)q'-\q - 1) if - (A, . . . , A- + 1, . . . , /. - 1, . . . , A), 

1 < fc < i, 

Uih + i)q''-\q - 1) if = (A, ■ • • , A - 1, • ■ • , A + 1, • ■ • , A), 

1 < fc < z, 



Proof: According to Property (v), the coefficients of the expansion of the product 
{e)Kf{e) into the basis {K^} are given by the intersection numbers of the scheme: 

(27) KFM)Kf{e) = Y^AjKhie). 

h 

Because the ordered metric is translation-invariant, we can assume that y = 0, so p^. ^ is 
the number of vectors z ^ j Q that satisfy z x for a fixed vector x 0. In other 
words, 

(28) z-x^ (0^ . . . , 0^ (ui, . . . , M„ 0, . . . , 0), 0^ . . . , 0"^) 

where the nonzero block is located in any of the n possible blocks, G F^, 1 < j < 
i,Ui ^ 0. 

(1) \h\ = I/I + 1. By the above we have that hj — fj for j ^ i and hi = fi + 1. 
Hence z can be chosen so that its A blocks of weight i annihilate the corresponding 
blocks of X, leaving one such block in any of the hi = A + 1 locations. Thus, 



The numbers pp. r are nonzero only in the three following cases. 



A + 1, /i = (A,---,A + i,...,A) 

otherwise. 



(2) |/i| I/I . The following numbers are easily verified by i 

h{q~2)q--\ h^f, 

(A + - i)q'-\ /i = (A, • • • , A + 1, • • • , A - 1, • • • , A), 1 < fc < 



^^•■^ 1 (A + i)(g-i)g'=-\ /i = (A,...,A-i,...,A + i,...,A),i<fc<^ 

otherwise. 

Other than these three cases, no other possibilities for h arise. 

(3) \h\ = |/| — 1. Now we should add to x one block of weight i in any of the n— I/I + 1 
all-zero blocks. Thus we obtain 

Pfu = (" - I/I + ^Wii - 1) h = (A, . . . , A - 1, . . . , A) 

and Pp, f — for all other h. 

To prove ( l25T l we now need to represent -P(e) as a linear combination of the KpiS. Using 
dlTT) we find that 

r 

|e|' = ^ze^ ^ (Jcrit^Ti. - ^LiKF^{e), 

i i=l 

hence 

r 

(29) P(e) = ^L,ifj.^(e). 
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The proof is now concluded by using dZTl l together with the intersection numbers computed 
above. I 

Along with the polynomials Ke below we use their normalized version Ke — KeJ -JvZ. 
The polynomials (Ke.^ e S A^.n) form an orthonormal basis of V . 

Denote by A,^, B^^, C^, the coefficient matrices of the normalized form of relation i 
The new matrix elements are given by 



A^[f, h] = + l)(n - K)qi-^[q - 1) if - (/i, ...,/, + 1, , 

h] = L,^{n -K+ l)M-^{q - 1) if/j- (A,. ..,/,-!,, 



J fr): 



(30) 



( UM-\q-2) 
q-l 



if h = f, 



L, 



B.[.f,h] = { 



q 



^{h + mq^' if/i = (/i,...,/fc + i,. ..,/,-!, 



1< fc < 



q 



•! fr)i 



1 <k <i. 

Let Vk C F be the set of polynomials of total degree < k. Let be the orthogonal 
projection of V onVi^. Define the operator 

Its matrix in the orthonormal basis has the form 



Bo 


Ao 


. . 


. 


Ci 


Bi 




. 





C2 


B2 .. 


. 







B. 



■ C, 

where the B^s are symmetric and Ci — Aj_^,i — 1 
and (|29] l. the matrix elements of are nonnegative. 
The matrix of S^, in the basis {Kf} has the property 

(31) VhS^[f,h]^VfS^[h,f] if,heAr,n). 



K. On account of property (v) 



(vii) (Explicit expression) 
(32) 



Kf{x) = q 



l/l'-l/l 



'i ; -^r— 2+1 ) 1 



i=l 



where fcj^ is a univariate Krawtchouk polynomial ( fT4l i. n,; — X]j=o ~ /i' 
and f.,xE Ar,„. This form of the polynomial Kf{x) was obtained in ifTOl (various other 
forms were found in [24iil4J). We remark that (|32] | can be proved by performing the Gram- 
Schmidt procedure (l22l l for monomials of higher degrees. It is known that the resulting 
system of polynomials is unique up to a constant factor once the polynomials of degrees 
and 1 together with the three-term relation jZSl l have been fixed, see 115, Theorem 3.4.9]. 

(viii) The matrix elements of the eigenvalue matrices of the association scheme H are 
given by P/e = Qfe = Kf{e). This follows from the previous property and (fT2l i because 
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the polynomials {Kf} form a unique orthogonal family on Ar,n with respect to the weight 

w{e). 

(ix) (Fourier transform representation) Let a; be a qth degree primitive root of unity, 
e,/ G Ar^n- Then 

(33) Kf{e)^ ^"'^ 



where shape(a;) = e. In fTol this relation is taken as a definition of the polynomials Kf. 
Under our approach, it follows from the well-known Fourier transform representation of 
the Krawtchouk polynomials ki{n, x) in the case r — 1 and Theorem B.ll 

(x) iChristojfel-Darboux). Let L C A,-.„ and define 

UL{a,e) ^ ^«7iX/(a)X/(e) (a,ee A,,„). 

Let S be the matrix of the operator S :Vn ^ Vn+i given by / ^ Pf, written in the basis 
{Kf}. The action of P{e) on Ul is described as follows: 

{Pie)-Pia))UL{a,e) = Y.'^f' E §[/, /i](X„(e)X/(a) - 

= E ^7' E m,hm,{e)Kf{a) - Kn{a)Kf{e)), 
feL heAr.r,\L 
the last equality justified by (ISTT l as follows: 

J2 vj'S[f,h]iK^,ie)Kf{a) - K^,ia)Kf{e)) 

f,h<£L 

= J2S[f,h]M{K4e)Kf{a) - i^/(e)if,(a)) 
=0. 

A particular case of the above is obtained when L — {f : \f \ < k}. The kernel Ul, 
denoted in this case by C/^, equals U^, = '^2 IKs(e)-^Ks(a), and we obtain 

(34) (P(e) - P{a))U,{a, e) - K,+i(e)^A^K«(a) - K^{ef A^K^+^{a) 
- ^ Q/(e)i^/(a)-i^/(e)Q/(a), 

/:|/I=K 

where Qf{e) = h]. This relation is called the Christoffel-Darboux 

formula. 

(xi) The generating function of the polynomials Kf is given by 

/ i=l j=l k=l 

In particular, 

^ Kf{e)=q^-6,,o- 

Remarks: 1. The polynomials J-iTe were considered in ll24l [T4l [TOl . However none of these 
papers constructed them from their definition as eigenvalues of the r-Hamming scheme 



16 



A. BARG AND P. PURKAYASTHA 



(to be more precise, Martin and Stinson 11241 mention this approach but pursue the path 
suggested in Theorem l3 . 1 I which makes explicit calculations difficult). Under the approach 
taken above, many properties of the polynomials K^. follow as special cases of the general 
combinatorial results of the previous section. 

2. Other generalizations of univariate Krawtchouk polynomials were considered earlier 
in ll37l [33)1 . These papers study biorthogonal polynomials for the weight given by the 
multinomial probability distribution, resulting in polynomial families different from the 
one considered above. 

3. Property (xi) implies a MacWilliams theorem for NRT codes. It was previously 
proved in [24. . 14 1 using different means. 

Theorem 4.2. (MacWilliams theorem in the NRT space). Let C d ~H and C-^ d H be 

two linear codes that satisfy 'YllLi ^iVi — for every x ^ C,y E . Let A{zq, z) = 
Tlie 111=0 weight enumerator of C and let A^(zo, z) be the same for . 

Then 

A-^{zq,Zi, . . . = |^yl(Mo,Ml, • ■ • 

where 

r j-i 
uq = zo + {q~l)^q'~'^Zi, Ur-j+i ^ zq + {q - 1) - q^^^z^, 1 < j < r. 

i=l k=l 

5. A Linear Programming Bound on Codes and OOAs 

In this section we prove one of our main results, an LP bound on the rate of codes. 
Let K = {Kf{e))f,eeA,^ „ be the eigenvalue matrix of the r-Hamming scheme, where for 
definiteness we are assuming the lexicographic order on the partitions. Let A — {A^), e S 
Ar,n, be a vector of nonnegative real variables, with the same ordering. Define two linear 
programs, 

(T) ■ V A ^ max ■ ^'^^0 ^ mill 

W ■ ^..^o^e-niax subject to 

(AkT>0 for all e (^^)^ = 1^"' ^ ^ 

and let LP(I) and LP(II) be their solutions. Let C{nr, M, d) be a code and C'{t, n, r, q) be 
an OOA of size M' . Then it follows from Delsarte's work ifTJl that 

M < LP{1), M' > iP(II). 

Replacing the LP problems by their duals, we obtain the following theorem. 

Theorem 5.1. Let F{x) = Fq + X^et^o ^eKe{x) be a polynomial that satisfies 

(35) Fo > 0, Fe>0(e^O); F{e) < for all e such that \e\' > d. 
Then any {nr, M, d) code satisfies 

(36) M<F{0)/Fa. 
Any OOA of strength t = d — 1 and size M' satisfies 

(37) M' > g"''Fo/F(0). 
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This result is essentially due to P. Delsarte. However, for the NRT space it was first 
stated by Martin ll24ll22ll and later rederived by Bierbrauer ifTol . The fact that the same 
polynomial gives a bound both on codes and orthogonal arrays is an easy consequence of 
Delsarte's theory, first mentioned in Levenshtein's work [|2T|. 

5- A. The bound. In this section we use Theorem lS.ll to derive a bound on ordered codes 
and arrays. Its proof uses a "spectral method" first employed in ID for the Grassmannian 
space and later used in |6 | to prove classical asymptotic bounds of coding theory. The gist 
of the method can be explained as follows. The polynomial F{e) is sought in the form 
F{e) = u(e)G^(e) where u{e) is a linear function that assures that F{e) < in (l35T l 
and G(e) is a function that maximizes the Fourier transform F{0). In the univariate case 
it turns out that a good choice for G is a delta-function at (or near) d. An approximation 
of the delta-function is given by the (Dirichlet) kernel which is its projection on V^. 
We therefore seek to modify the operator so that becomes its eigenfunction with 
eigenvalue 6'„, express the bound of Theorem 15. II as a function of 9k, and optimize on k 
within the limits ( |35] |. The reader is referred to Q for a more detailed discussion of these 
ideas. 

Theorem 5.2. Let k be any degree such that P{e) < X^-ifor all shapes e with \e\' > d, 
where Xi is the maximum eigenvalue of Si and d > 1 is an integer 
Let C G H be an [nr, M, d) code. Then 

(38) ^^^ Ar5.An-.W-lT fn-^ 

Scvurn - Ak V'*. 
Let C be a (t ~ d — l,n,r,q) OOA of size M. Then 

(39) M>'^"' iS.r,rn^K) 



□ 4r5erit(n-«)(g'' 

Proof : Consider the operator that equals Sk, on Vk_i and acts on a function ip e 

14\K-i by 

f--\f\=i^ 

where e/ > are some constants indexed by the partitions of weight k (their values will 
be chosen later). The matrix of in the orthonormal basis equals 



E 

where E = diag(£/, |/| = k) is a matrix of order {'^^.L^^) ■ Let mbe such that + to/ > 
0. By Perron-Frobenius ifTTl p. 80], the spectral radius p(Tk + ml) is well defined and is 
an eigenvalue of (algebraic and geometric) multiplicity one of + ml. Moreover, again 
using Perron-Frobenius, 

p{Sk-i + ml) < p{Tk + ml) < p{Sk, + ml). 

Then 

(40) Xk-1 <0n<XK 

where 9k. = p{Tk)- Let G > be the eigenfunction of Tk with eigenvalue 9k. Let us write 
out the product P{e)G in the orthonormal basis: 

F(e)G = 5„G + G,A„K„+i=^?«G+ ^ SfGfKf + GkAkKk+l 

f--\f\=i^ 
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where Gk, is a projection of the vector G on the space Vk\Vk-i. This implies the equality 

i:\n=.Gf{efkf + Qf) 



G 



P{e) - 



where Qf{x) is defined after (|34] |. Now take F{x) — (P{x) — 0^)G'^{x). Let us verify 
(|35] |. Since multiplication by a function is a self-adjoint operator, we obtain 

Fo = (^^, 1) = ( ^ GfiefKf + Qf),G) = ^ G)ef > 0. 

Using (l24l l one can easily check that Fe>Q for all e 7^ 0. The assumption of the theorem 
together with (gOll implies that F{f) < for |/|' > d. Hence 



{j:\fi=.Gf{efKf{0) + QfiO)) 
1 ^ {efKf{0) + Qf{Oj) 



E 



where in the last step we used the Cauchy-Schwarz inequality. Computing the minimum 
on £/ , we obtain 

(41) M^j^(^^EQf(^)^- 

Next, 

\f\=K /:|/|=K h:\h\ = K+l 

Let h — (/i, . . . , /i + 1, . . . , fr) for some i,l < i < r. Then using ([1) we find 



A, [/, /i] ^ v/(/» + l)(n - k)<7*-i(9 - 



Thus we have 



1 

E Q/(o)V^= E E("-^)(i-7^)^/ 

|/|=K |/|=K« = 1 



= (n - K)r(5ci-it ^ w/ = (n - K)r(5ci-it f ](<?''- 

Substitution of this expression into ( l4TT l concludes the proof of ( l38T l. The bound (|39] l 
follows bv (|37]|. I 
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5-B. Spectral radius of Sr. In this section we derive an asymptotic lower bound on the 
spectral radius of S^- This estimate will be later used to optimize the bound dSSl) on the 
choice of k. 



Theorem 5.3. 



lim — > max A(ti , . . . , Tr), 



where 



A(ti, . . . , Tr) =Y.L, (2v/(l - rMq - 



(42) +^q-2)n{q^~q^-') + 2^-^^Y.^' 



k=i 



TkTiQ 



i+k 



To prove this theorem, we will bound below the largest eigenvalue A^; of the matrix S^- 
For any real vector y we have 

{y,y) ' 

We will construct a suitable (0, l)-vector y. Its coordinates are indexed by the partitions 
arranged in the increasing order of their length /i and lexicographically within a block 
of coordinates for each value of /i, < /i < n. Let y ~ {yo,yii ■ ■ ■ lUn) where = 
(2//,I/I=m)- 

Let /, I/I = /i, be a shape vector. For an integer J consider the set 

r 

= ^1^{JJ) - {(/l + K + jl,. .. ,fr+ ir) ■■^jt^Q\\ji\< J,i = l,-- - .r} 

4=1 

and denote m = \J-fj,\. Next, let 

iy^)h = i{he T^) 

for /i = K + l — s,...,K where s will be chosen later, and y^ — otherwise. 

In the next two lemmas we derive a lower bound on the part of the product y'^S^y that 
involves only the rows of that correspond to the shapes / of length /i. Let 

Eh = {{hi, . . . , /ifc ± 1, . . . , /ij T 1, • • • , K), l<k<l<r} 

be the index set of the nonzero off-diagonal elements in the row in which is indexed by 

h = {hi, . . .,hr). 

Lemma 5.4. Lete = argmin^^^^ (X;ge£;,u{h} an-^^ef V'^ ^ Eg-S^[e,5]- Then 

y^Bf^yf, > '(/'^m(l - o,„(l)). 

Proof: I. Since \h\ = fj, for every h e J-fj,, the quantity IJF^I equals the number of ordered 
partitions of into at most r parts, each part bounded between — J and J, or the number 
of ordered partitions 



Jr = J2J^, 0<j,<2J,i^l, 
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The number of such partitions is given by LIT, p. 1037]: 

Writing this expression as a polynomial in J, we find the coefficient of J*"^^ to be 

^|(-i,-(;:),.-.r- 

Since this is always positiv^l we conclude that m is a degree- (r — 1) polynomial in J; in 
particular, if J ^ oo, then also m ^ oo. 
II. We have 

To bound yJ^B^y^j^ below we estimate the difference between the above sum and the sum 
of all the nonzero elements of in the rows /i G JF^ which is obtained by replacing the 
range of column indices g above with g G U {/i}. Therefore, for a given h G let us 
estimate the number \Eh\T^\ of nonzero entries in B^[h, ■] not included in the sum. Let 
/ = (fi, . . . , fr) and let h be of the form /i =(..., + J, ...) G JF^i for some 1 < k < r. 
Consider the column indices g ^ Eh given by 

(43) g = (/i + At - K + ji, . . . , /fc + J + 1, . . . , /, + J, - 1, . . . , + j^) 

for any k ^ I ^ {1, . . . ,r} . For any pair ft,, g of this form, B^[h, g] ^ but g ^ T^^. The 
number of shapes h that result in shapes g of the form (|43] | equals the number of ordered 
partitions of — J into at most r — 1 parts of magnitude < J; equivalently, this is the number 
of ordered partitions 

J(r - 2) = j2 + • • ■ + Jr, < 2, < 2 J, i = 2, . . . , r, 

which equals n_(- = 7r(r — 1, 2 J, J(r — 2)). 

Next consider the row indices h= (. . . , /fc — J, . . . ) G J-^, and column indices g ^ Eh 
given by 

9 ^ ifl + fJ- - K + ji, . . . Jk - J - 1, . . . , fl + jl + 1, . . . , fr + jr) 

which again account for B^[h, g] ^ and g ^ T^^. The number of such shapes h equals 
the number of ordered partitions of Jr into r — 1 or fewer parts < < 2 J. Denote 
this number by n_ = 7r(r — 1, 2 J, Jr). Note that as J ^ oo, both !!-)_ and n_ grow 
proportionally to J''"'^. 

It is easy to verify that Eh \J-f^ 7^ if and only if h and g are of the described form. 
Observe that by ( |30] |. \Ef \ = — r. We then obtain 

5] > V'mC"! - ri^^ - 0(n+ + n_)) = ^^m(l - o„,(l)). 
The lemma is proved. I 

' To prove positivity, observe that the numbers Sr,m = '^\^q {\){''" ~ Si)™ satisfy the recuiTence 

Sr,m = r'^Sr,m-2 + 4r(r - l)5r-2,m-2, 3 < m < r — 1 
and then use induction to prove that 5,-.m > (< 0) according as r — m = 1 or 3 mod 4. 
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We now consider the part of the product y^S^y that involves the matrix C^^, fj, = k — 
s + 2, . . . , K. For a shape h let 

Dh = {{hi,...,hk-l,...,hr),l<k<r}. 

The proof of the next lemma is very similar to the above proof and will therefore be omit- 
ted. 

Lemma 5.5. Let e = argmin/^^^^ (X;^^^,^ G^[h,g]) and let (j)^, = 

ylCfiy^i-i > ^^"^(1 - Om(l))- 
To complete the proof of Theorem l5.3l compute 

1 J. 
Ak > — y S^y 
ms 



= — ( E ylB,y, + 2 ylc,y,~i 

^ — K-\-l — S fl — K-\-2 — S 

>-( E ^^. + 2 ^ -o„,(l)) 



fl—K-\-l — S /.i — K + 2 — S 

>^* + 2^r(l-o„(l)), 
s 

where ip* ((j)*) is the smallest of the numbers (0^) above. Note that both ■(/;* and 0* are 
nonzero. Now let n — > oo, k = rn, and let us choose / in the definition of JF^ to be of the 
form / = (/i, . . . , fr), fi = riTi, 1 < i < r. We assume that none of the r/s approach 
as n grows. Take s — o{n), s oo. Using (l30l) . and letting J = 0(71), J — > 00 we get 



lim ^ = V - 2)r,(g'- - 9^-1 ) + 2^^-^ V V^^) 

i=l ^ fc=l 

lim — = V [I - T)n{q - l)q^-K 



n — *oo Tl 

i=l 



Then since n/n t, 



lim — > lim — — > A(ti, . . . ,Tj.). 

n— too n n— too mSU 



The theorem is proved. 

5-C. Asymptotic estimate for codes and OOAs. Theorems 15 .21 and 1531 together enable 
us to prove one of the main results of the paper. 

Theorem 5.6. Let Rij>{d) be the function defined parametrically by the relations 

(44) R^r)^-(h,{T)+T\og^^—^ 

r \ ^ q — 1 

(45) 6{t) = S„n ~ - max A(ti, . . . , t^), < t < 1. 

r Ti>o 

Then the asymptotic rate of any code family of relative distance S satisfies R < R]^p{d) 
and the rate of any family of OOAs of relative strength 6 satisfies R> 1 ~ i?Lp(<5). 
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To prove this theorem, take the logarithms in ( I38l l and pass to the Umit as n oo. 
Using the standard asymptotics for the binomial coefficient, we find that the code rate is 
bounded above by the right-hand side of ( l44b . The condition P(e) < A„_i of the Theorem 
15. 21 will be satisfied for large n if 

Sail - S < 

rn 

This defines the function in ( |45l l. Thus, the proof is complete. 

Remark. For r = 1 this bound reduces to the linear programming bound on the rate of 
codes in 12. Just as that result, the bound of this theorem improves upon the asymptotic 
Plotkin bound for large values of the code distance. 

6. The Case r = 2 

In this section we prove a bound for codes in Q^ " which improves upon the general 
result of the previous section. The improvement is due to the fact that in the case r = 2 it 
is possible to work with the polynomials Kf{e) in their explicit form, and base the bound 
on the behavior of their zeros instead of the spectral radius of the operator S'^ . Namely, let 
/ = (/i, /2), e = (ei, 62). From ^ we have 

Kf{e) = q^^k^^ (n - 62, ei)fc^^ (n - /2, 62). 

We also have 

P(e) ^n(2- ^-tij - ei - 2e2. 

We will use the following properties of the polynomials kg whose proofs are found for 
instance in |21|. Let Xi(n, s),i — 1, . . . ,she the roots of kg in the ascending order. Then 

(46) Xi(n — 1, s) < Xi{n, s) < Xi{n — 1, s — 1) < Xi{n, s — 1) < Xi^i(n^ s), 

1 < s < n, i = 1, . . . , s — 1. 

Let n — > 00, s/n y. Then 

(47) hm = 7 y = y - - V 9 - l)y 1 - 2/ • 

Ti-^oo n 1 1 Q 

The Krawtchouk polynomials satisfy the recurrence 

(48) ks{n, x) — kg{n — 1, x) + (9 — \)kg_i{n — 1, x) 
and a Christoffel-Darboux formula of the form (|34] | 

(49) q{x ~ y) '^'^k^^O)^^ = 'j^ikh+i{y)kh{x) - kh+i{x)kh{y)). 

Remark: Properties (l46ll-(l49ll are usually stated for integer n. This is related to the 
fact that the polynomials ks{n, x) represent the eigenvalues of the Hamming association 
scheme. As pointed out to us by M. Aaltonen |3|, it is possible to prove these properties 
for any n G relying on the generating function of the Krawtchouk polynomials. 

The main result of this section is given in the following theorem. 

Theorem 6.1. The asymptotic rate of any family of codes of relative distance S satisfies 
R < <f>((5), where 

$(5) = min V2{t2 + /i,(ti) + (1 - ri)/i,(-^)}, 
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where the minimum is taken over all ti , T2 that satisfy 

< ri < ((7 - l)/q\ 0<T2<{q^ l)/q, 
7(t2) + (2 - 7(t-2))(1 - T2)7(ti) < 25. 

The asymptotic rate of any family of OOAs of relative strength 5 satisfies R > 1 — ^{^)- 

The remainder of the section is devoted to the proof of this result. We note that the 
polynomials Kf{e) are formed as products of two Krawtchouk polynomials. A similar 
situation arose in [l] which dealt with the Johnson association scheme whose (second) 
eigenvalues are equal to a product of a Krawtchouk and a Hahn polynomial. Therefore, we 
adopt some elements of the analysis in [T| in our proof below. 

In quest of an LP bound, we require a polynomial F{e) = F{ei,e2) that satisfies 
conditions (l35l l. Consider the polynomial of the form 

(50) F{e) = (P(e) - P{a))iUL{a, e)f 

for some a — {a, (3) and a subset L. For brevity below we write Sfh instead of Sk[/, h] 
and denote L — A2.„\i. We find 

Fo = {F, 1) = ((F(e) - P(a))t/L(a, e), U^a, e)) 
feL f heL 

In order to ensure that Fq > we will choose L and a so that 

(51) Kh{a)<{) if heL; Kf{a)>0 if/sL. 

Let s = (si — 1, S2) G A2,„ be a shape that satisfies {(si — 1, S2 + 1), (si, S2)} C A2,„. 
Let a = (a, (3) satisfy 

(52) (3 = xi{n — S2, si), xi{n — (3, S2 + I) < a < xi{n — f3, S2), a + 2/3 < d. 
For any /2, < /2 < S2 + 1 denote by 0(/2) the degree such that 

(53) xi{n - f2,<f>{f2) + I) < 13 < xi{n ~ f2,4>{f2)). 
By ( |46] |. (/){■) is well defined and implies the following: 

[{xi{n - u,w) > P) ^ [w < (t){u))], [{xi{n~ u,w) < (3) ^ (w > + 1)]. 

We choose the region L to be given by 

i=(/i,/2:/2 = O,...,S2;/i=O,...,0(/2)). 

For the moment this choice is not unique because there are many possibilities for s. This 
ambiguity will be later removed by optimizing the bound on the choice of s. 

To argue about the sign of Fq we need to establish some properties of the region L. 
First, we claim that for a fixed /2, 

(54) ^{f2) ~ 1 < Hf2 + 1) < 0(/2). 

Indeed, by ( |46] l. 

P < a;i(n- /2X/2)) < xi{n - f2 - l,9i(/2) - 1) 
which implies the left-hand side of ( |54l ). On the other hand, 

I3>xiin- /2, (/-(/a) + 1) > - /2 - 1, <^(/2) + 1) 
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which implies the right-hand side. 

The values of /, h for which S/h ^ are given in dSOl) . In particular, if / G i, then 
the set H of the shape vectors h that index the nonzero matrix elements of S and that lie 
outside the region L is as follows: 

i^ = {(0(/2) + 1, h). /2 = 0, 1, . . . , 52} U {(/i, S2 + 1), /i = 0, 1, . . . , Si - 1}. 

The region L and the corresponding set H are shown in Fig. [T| By our choice of the 
parameters, 

/c/,(n-/3,a) >0 (0 < < sa), 
fcs2+i(»^-/3, a) < 0, 
kf,{n ~ /2,/3) > (0 < A < 0(/2), < /2 < S2 + 1), 

fc*(/2) + l(«-/2,/3) <0 (0</2<52 + l). 



m 



All points X belong to the set H 
All points • belong to the boundary of L 

L 




Figure 1 . The region L 



Then 

(55) K(^f^j,_){a)^qf^kf^{n- (},a)kf,{n- f2,f})>Q, f e L, 

(56) K(^(^f^)+ij^){a) ^ q^^kf^{n- /3,a)k^(f^)+i{n~ f2,f3) <0, < h < S2, 

(57) (a) = q''+^ks,+i (n - /?, a)fc/, (n - S2 - 1, /3) < 0, < /i < si - 1. 
Thus, 

(58) Kf{a) < (/ e L), Kf{a) > (/ G L). 
This proves that Fq > 0. 
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Let us show that Fe > for all e. For this rewrite F as follows: 

F{e) ^ {{P{e) - P{a))UL{a,ef 

= E f E Sjh{Ku{e)KM - K,{a)Kj[e)) ^ 

= ( E E "^-^^^-^^"^ - J2 ^/(^) E ^''^"^'^^^' ) ^ -^3(«)-^g(e) 

- E ^i^/(e)i^.(e)E^^^- 

By ( |24] |. the products KhKg and KfKg are expanded in the basis {A'/} with nonnegative 
coefficients. Moreover, the other terms in the above formula also have the needed signs on 
account of (|55]|-(|57]|. This establishes our claim. 

Finally, because of the third condition in (|52] |. F{e) < for all e with |e|' > d. 

We are now able to formulate the bound on codes and OOAs. 

Theorem 6.2. Let C be an {2n, M, d) code C C H{q, n, 2). Then 

(59) M < r-^-r Vs, 

where s — {si — 1,82) satisfies {(si — 1,S2 + l),(si,S2)} C ^2,n and a — (a,/?) is 
chosen tofiilfill conditions (l52ll. 

Let C be a (t = d — 1, n, 2, g) (90A of size M. Then 

(60) M > — _ ^ _ _ _ + i)2(^ _ i)3(c, + 2/3) ■ 

Proof: Let us compute Fq = (F, 1). Denote cti = (si — 1, S2 + 1), o'2 = (si — 2, S2 + !)• 
By dSSl ) and ( l26b we have 



/6i /iGL 

> - •^"^"'^ (5^,cri ATo-i (a) + S's,o-2Aro-2(a) + S'5,(si,s2)-f'^(si,s2) (a)) 



Ks{a) 



q2^{('^ + 1)^-1 («) + + 1)(9 - 1)^.2(0) + si{q + l)A-(,^,,,)(a)) 

(61) - _ ^.(«)(^2 + l)g-- + ^ ^^^^^^^ _ _ ,^^^), 

Let us now evaluate UL{a, 0). 

S2 <t>{f2) 

UL{a, 0) = E = E 9^'^^^ (" E % (" - /2, 

/ei /2=o /i=o 

Let us bound above the last sum. We shall prove that 

0(/2) 0(^2) 

(62) (" - /2, /?) < q'^-^' -S2,P). 

/i=0 /i=0 
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Indeed, using (I48b . we obtain 

0(/2) Hf2) 

E % (" - /2' = E (^/i (" - /2 - 1. /3) + (9 - l)%-l(^ - /2 - 1, P)) 
fi=0 /i=0 

0(/2)-l 

= - /2 - 1,^) + g E /2 - I./?)- 

/l=0 

Recall that (/)(/2 + 1) = (^(/2) or </i(/2 + 1) = 4'if2) - 1- In the former case, 

0(/2)-l 0(/2 + l) 

/i=0 /i=0 

in the latter, k^(f2){n - /2 - l,/3) = fc0(/2+l)+l(?^ - /2 - 1, /3) < on account of 
and (l46T l. Repeating this procedure S2 — /2 times, we arrive at (|62] |. 
Note that 0(s2) = si — 1. Therefore 

^^2 Si- 1 

UUa,0)<q'' J2 kf,in-(3,a) ^ % (n - S2 , /?) . 

/2=0 /l=0 

By (|49ll, and (|23]l we have 

Si(T')(9-l) 
> kf,{n-S2,P) = afn-S2\ ks,^i{n - S2, P) , 



(s2 + I)(fc,s2+i(n',0)fcs2(n',a) - fcs2+i(n',a)fc52(n',0)) 



gQfcs2(n',0) 

/2=U 

qa 

where n' ^ n — (3 and W{x) = kg^+iin', x)/ks2 {n' , x). Using these expressions, we can 
bound UL{a, 0) as 



{S2 + l)(n - S2 ~ Si + l){q- 1) 
q'^a/] 

Hence using (l36l) . dSOl) . and ( |6T1 ) we can write 

(S2 + l)(n - S2 - Si + - + 2/3) (W^(0) - T4^(a))2 



Since W{0) = (g — l)(n' — S2)/(s2 + 1) > > 1^(0;) > —00 as a ranges in between the 
bounds in ( |52] |. it is possible to find a such that VF(a) = —W{0). With this choice and (HJ 
the last expression turns into ( l59l ). The estimate ( l60l ) follows from djTl l. I 

The proof of Theorem 16. H is obtained by passing to asymptotics in ( [59] ). Namely, let 
n —^ 00, d/nr —^ d, si/n — > ti, 52/?^ — > T2. By ( |47| l. 

P a 
limsup- = 7(n)(l - T-2), limsup - = 7('r2)(l ~ 7(ti)(1 - T2)). 
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r=2 q=2 




1=2 q=3 



0.3 0.4 0.5 0.6 0.7 




Figure 2. Bounds for codes. 



Computing the logarithm on the right-hand side of i5% . we observe that the only term of 
exponential growth arises from Vs . Using standard estimates we obtain 



log„ Vs = log 



S^ -1 



n — si 

S2 



+ S2-1 



< n{T2 + hgin) + (1 - Tl)hg{j^)} 



The tightest bound is obtained by computing the minimum of this expression on ti,T2. The 
range of the variables ri , T2 is obtained on observing that n{q~l) jq^ and 77.(17—1) jq are the 
maximizing values of si, S2 for large n (by a direct calculation from the above expression; 
or, specializing from a general result in |34|). The third restriction in the statement of the 
theorem is implied by a + 2/3 < d. This completes the proof. 

Asymptotic bounds for ordered codes are shown in several plots in Fig.|2] 

Remark. The NRT metric is an example of a wide class of metrics on the set termed 
poset metrics following Brualdi et al. 1 1 2 1 . To define a poset metric, consider a partial order 
on the set [1,2,..., N\. An ideal in the order is a subset closed under the ^ relation. 
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The -weight of a vector x e Q is the size of the smallest ideal that contains the 
nonzero entries of x. For instance, for the NRT weight, the relation -< can be defined as 
-< (*2, ^2) iff ii = 12, ji < j2, where ii,i2 are the indices of the block. A dual 
order Py is formed of the same set of chains as the order but with the signs reversed 
within each chain. In ll24l [35l IT4l . and in our paper, the poset duality is realized as the 
C C H , C H convention. One of the main questions that arises in this context is to 
characterise the association scheme that arises from the order and in particular, to derive 
the Mac Williams-type relations. Partial orders that give rise to a univariate Mac Williams 
relation have been described by Kim and Oh 1 19 1. On the other hand, rather little is known 
about the multivariate case of which the NRT space is an instance. 

Acknowledgment: A.B. is grateful to WilHam Martin for calUng his attention to the 
problem of code bounds for the NRT space. The authors are also grateful to M. Aaltonen 
for a useful discussion of his work HI. 
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